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Abstract 



We study the dynamical and statistical properties of turbulent cross-helicity 
(correlation of the aligned fluctuating velocity and magnetic field components). 
We derive an equation governing generation and evolution of the turbulent cross- 
helicity and discuss its meaning for the dynamo. Using symmetry properties of 
the problem we suggest a general expression for the turbulent cross-helicity. Ef- 
fects of the density stratification, large-scale magnetic fields, differential rotation 
and turbulent convection are taken into account. We investigate the relative con- 
tribution of these effects to the cross-helicity evolution for two kinds of dynamo 
models of the solar cycle: a distributed mean-field model and a flux-transport 
dynamo model. We show that the contribution from the density stratification 
follows the evolution of the radial magnetic field, while large-scale electric cur- 
rents produce a more complicated pattern of the cross-helicity of the comparable 
magnitude. The pattern of the cross-helicity evolution strongly depends on de- 
tails of the dynamo mechanism. Thus, we anticipate that direct observations of 
the cross-helicity on the Sun may serve for the diagnostic purpose of the solar 
dynamo process. 

Subject headings: Sun: dynamo — Sun: interior — Sun: magnetic topology — 
turbulence 
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Introduction 



The cross-helicity conservation law has been estabhshed by |Woltjer (1958) and Moffat 



(1969). It states that in perfectly conductive media the parallel components of magnetic 



and velocity fields do not interact. More precisely, the conservation law states that if the 
velocity field of plasma, U, and the magnetic induction vector field, B, are confined in the 
volume V then the integral / = /y U • BdV is time-invariant if the dissipative processes are 
absent. 

The "Alfvenic" state of magnetohydrodynamical (MHD) turbulence, in which the 
fluctuating magnetic and velocity fields are aligned, is commonly considered as a preferred 



state in MHD relaxation processes, (Woltjer, 1958; Matthaeus et al., 2008). Such state is 



also accompanied by quenching of the turbulent magnetic field generation effects (e.g., the 
a-effect) due to the back reaction of large-scale magnetic fields on the turbulent motions. 
Therefore, the cross-helicity can be a useful quantity for diagnostics of non-linear turbulent 



dynamo processes in the solar convection zone, (see, e.g., Kleeorin et al., 2003). Yoshizawa 



(1990) and Yokoi (1999) considered the cross-helicity as a part of the dynamo mechanism 



in turbulent astrophysical fiows. The role of the cross-helicity conservation for this type 



of dynamo was explored recently by Sur & Brandenburg (2009). The results of Riidiger 



et al. (2010) suggested that the cross-helicity parameter that could be observed on the solar 



surface is a manifestation of solely near-surface physical processes. On the other hand. 



the results of numerical simulations (Mason et al., 2006; Boldyrev et al., 2009) showed 



the existence of rapid local alignments of turbulent velocity and magnetic field in the 
presence of spatially nonuniform pressure and kinetic energy of turbulent fiuct nations. It is 
notable that these locally aligned turbulent velocity and magnetic field structures spatially 
dominate even in the case when the mean cross-helicity is zero. Some initial attempts to 



measure the cross-helicity from SOHO/MDI data ( [Scherrer et al. , 1995|) were made by [Zhao 



eraD(l20TTl). 



In observations of the cross-helicity at the solar surface the processes of the "local" 
and "global" alignments may be mixed. To isolate the effect due to the local alignment we 
have to choose appropriate spatial and temporal scales for averaging. In this paper, we 
show that the large-scale distribution of the averaged cross-helicity can be affected by the 
solar dynamo operating in the deep convection zone, although the scale separation problem 
remains a difficult issue of the theory. 



The theoretical approach to study the solar dynamo is based on the mean-field 



magnetohydrodynamics (Krause & Radler, 1980). In the following part of the paper we 



propose a general expression for the turbulent cross-helicity by the use of the transformation 
symmetry properties. Then we estimate the coefficients in this expression by means of 



the T- approximation in the turbulence theory (Orszag, 1970; Vainshtein & Kitchatinov 



T9831 IBlackman fc Brandenburg[ [20021 [Field fc Blackmanj [2002[ [Blackman fc Fieldj [2002 



Brandenburg & Subramanian, 2005). Finally, we illustrate the major contributions to the 



cross-helicity evolution using two dynamo models of the solar magnetic cycle. 



2. Transformation symmetry properties and cross-helicity 



Following Krause & Radler (1980), we split the physical quantities of the turbulent 



conducting fluid into the mean and randomly fluctuating parts. The mean parts are defined 
as ensemble averages of the corresponding properties. Magnetic field B and velocity U are 
decomposed as follows: 

B = B + b, U = U + u. (1) 

Hereafter, we use capital letters with a bar above for the mean-field (large-scale) properties 



and the small letters for the fluctuating (small-scale) parts. As shown by Krause & Radler 
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(1980) and Raedler (1980) the general structure of the mean electromotive force of turbulent 
flows, ^ = u X b, can be reconstructed by using the transformation symmetry properties of 
the basic physical quantities and their products involved in the problem. We deal with the 



cross-helicity pseudo-scalar u • b in the similar way. 

Any pseudo-scalar can be expressed either as a tensor product of a tensor with 
pseudo-tensor, or as a scalar product of vector and pseudo- vector. General expression 



for u • b may have a fairly complicated structure. We restrict ourselves in studying the 
effects that can be important for the solar convection zone dynamics and solar dynamo. 
Therefore, we have to take into account the stratification of the thermodynamic and 
turbulent parameters. Their contribution to the solar dynamics is given by gradients 
of the turbulent kinetic fiuctuations and magnetic energy density, e.g., V (pu^l and 



b^ 

V [ - — I (all quantities are ordinary vectors). Note, that the mean density is suflicient 
to describe the effects of the thermodynamic stratification in the polytropic atmosphere. 
This assumption will be used hereafter in the paper. Also, we have to take into account 
effects of the global rotation (pseudo- vector Ct) and large-scale shear, ViUj (tensor). This 
tensor can be decomposed as follows: V^t/j = \sijpWp + (VU)^^^.-^, where, W = (V x U), 
(pseudo- vector) and a strain tensor, (VU)^^^.-^ = | {^iUj + VjUi). We assume that 
u • b = in the absence of large-scale magnetic field and take into account effects of the 
large-scale magnetic field, B (pseudo- vector) , and its spatial derivatives, ViBj. Similar to 
the large-scale shear fiows we decompose the contribution of the non-uniform magnetic 
field into effects of large-scale electric current J = V x B//xo (vector) and magnetic strain 



tensor, (VB)^.^.-^ = | (V^Sj + VjS^) ( pseudo-tensor). By analogy with Raedler (1980), 
and also assuming the scale separation {£^ Tc} <C {!/, T} of the typical small-scale (^, Tc) and 
the large-scale (L, T) spatial and temporal variations of the velocity and magnetic fields, we 
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write a combination of above effects as follows: 



u.b = |(B.V)(^«.pu^ + «.^j (2) 

where, /^i_5 are coefficients that needs to be determined). Having in mind the 
transformation symmetry properties of the quantities that are involved in the problem, we. 



b2 



of course, can construct the pseudo-scalar u • b in many other ways. In principle, the Eq. 
([2]) may include the terms like, G • 12, U • 12 (and similar), or the terms like, U • B, B • J 
(and similar), or even higher order combination of the physical quantities. Later, we will see 
that such terms do not appear in our analysis of the momentum and induction equations of 
the fluctuating velocity and magnetic flelds if we restrict ourselves to the case of the weak 
large-scale magnetic fleld B. In this sense the Equation ^ may be incomplete and should 
be considered with cautions in analysis of the observational data. 



The physical interpretation of the first term in Eq. ^ was given by Riidiger et al. 



(2010). Consider the turbulent medium permeated by the large-scale field B. Convective 
elements rising, with velocity u expand, (V • u) > 0, and induce a fiuctuating magnetic 
field, b ^ —TcB (V • u) (sign is anti-correlated with sign of (V • u)). The same is valid for 



descending and contracting convective elements. In the anelastic approximation (Gough 



1969), V • pu = and (V • u) = — (u • V log p). Therefore the sign of u • b is opposite to the 
sign of B ( because logp < 0). The other terms in Eq. ^ are less easy to interpret. The 
effects of density fiuctuations are not taken into account in this consideration. For weakly 



U2 



compressible (subsonic) convective fiows, <C 1 {Cg is the sound speed), the contribution 



r^2 



U2 



of the buoyancy forces to the cross-helicity is proportional to ^ i^ey^^c (g * B), where g is 
the gravity acceleration. This effect may be of the same order magnitude as the first term in 



U2 



Eq. (|2|), particularly close to the surface where can be significant. However, we do not 



CI 
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consider the density fluctuations effects, and limit our study to the anelastic approximation, 
which aUows us to investigate the turbulent cross-helicity analytically. In what follows 
we estimate the one-point cross-helicity correlations using two different approaches in the 
mean-fleld theory, based on the so-called r- approximation of the turbulence theory, in 
which third-order correlations are approximated by a dissipative term with a characteristic 



dissipation time Tc (Orszag, 1970; Vainshtein & Kitchatinov, 1983; Kleeorin et al., 1996 



Rogachevskii fc Kleeorin[ |2000[ [Blackman fc Field] |2002[ [Blackman fc Brandenburg[ |2002 



Brandenburg & Subramanian, 2005). Using the r- approximation we deal with small 



deviations in the turbulent state resulting from large-scale magnetic flelds and flows. In 
this approximation it is assumed that the background turbulent state (in the absence 
of the large-scale flelds) has a priori deflned statistical properties. Our flrst approach is 
based on the cross-helicity evolution equation and averaging turbulent properties in the 
physical coordinate space. It is simpler but less accurate than the second approach that 
uses a two-scale Fourier transformation for equations governing the evolution of fluctuating 
velocity and magnetic flelds. 



3. Mean-field theory of the cross-heUcity 



3.1. Calculation of u • b by averaging in the physical coordinate space 



As pointed out by |Yoshizawa et al. (2000) the evolution equation for the cross-helicity 
is useful for exploring mean-fleld dynamo properties based on the cross-helicity effects (also 



see, 



Sur & Brandenburg, 2009). We start with the ideal MHD induction and momentum 



equations written in a coordinate frame rotating with a constant angular velocity Q. : 



~dt 

pdtl^ + p (U- V) U + 2p (n X U) 



V X (U X B) , 

(V X B) X B 

Mo 



(3) 
(4) 
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Here, U and B are the velocity and magnetic field vectors, p is the gas pressure, "i/g is the 
gravitational potential and T is the external random force. Calculating the scalar products 
of the first equation with U, and the second equation with B/p, and summing them up, we 
arrive to: 

ai(U-B) + 20-(UxB) = -V-|u(U-B) + B(^*, + ^-^-H!)| + (B.f')(5) 

\/p f 1 P 

Here, we assumed the polytropic pressure law, p ^ (this gives — = V 



Eq.([5]) describes the cross-helicity conservation of barotrophic fiuids (Woltjer, 1958) in the 
inertial frame (case = ) and for a potential external force, T = —Vcj). The conservation 
law states that dt /y U • BdV=0 in a volume V if the normal components of U and B 
vanish at the boundary. In the turbulent media the cross-helicity is not conserved because 
nothing prevents it to cascade from large scales to small scales where it is dissipated. This 



problem was recently discussed by Riidiger et al. (2010). 



As discussed above, in our study we neglect the density fiuctuations effects and adopt 



the anelastic approximation, divpC/ = divpi^ = 0, (Gough, 1969). Substituting Eqs.rtl^ in 



Eqs.([3]) and Eq.Q and averaging them over an ensemble of fiuctuations we get 

dtB = V X (5 + U X B) , ^ = loHb (6) 

— / — N (V X B) X B , _ 

pdtU + 2p{nxU) = ^ ^ p(U-V)U-pV^^ (7) 

Po 

(V X b) X b 



- Vp+ ^ p(u-V)u. 

Po 



Similarly we derive the following evolution equation for u-b after subtracting the equation 
for the mean-field cross-helicity U • B from Eq.(^5h: 



dt{u'h) = -V-J'^-^-(2n+(VxU)) + ^V,r,, -^^^ (8) 

P 

1 / 

Tij = puiUj bibj - 5ij— , (9) 

Po y ^ J 

= TJu^ + BilMU-B—. (10) 



Bi (u-b) 
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Here, following the r-approximation approach (Vainshtein & Kitchatinov, 1983; Kleeorin 



eFalj p96l [Blackman fc Fieldj [20021 [Blackman fc Brandenburg[ [20021 in the Eq.Q, we 



replace the third-order correlations of the fluctuating parameters, and the terms b • f 



with a relaxation term 



(u-b) 



, here Tc is a typical time scale of the turbulent motions. 



Approximating these complicated contributions with the simple term 



(^) 



has to be 



considered as a questionable assumption. It involves additional assumptions (see Radler & 



Rheinhardt, 2007), e.g., it is assumed that the second-order correlations in Eq.(8) do not 
vary signiflcantly on the time scale of Tc. This assumption is consistent with scale separation 
between the mean and fluctuating quantities in the mean-fleld magnetohydrodynamics. 
The reader can flnd a comprehensive discussion of the r -approximation in the above cited 
papers. 

(J 

In our study we discard the cross-helicity flux in the further analysis. This can be 
partly justifled if we average over a sufliciently large volume of the fluid in which the vector 
fleld T is conflned. For example, if the averaging is carried out over an ensemble of the 
turbulent flelds, then we assume that this ensemble is sufliciently representative, and that 
all members of the ensemble are conflned and uniformly distributed in the volume. 

In addition, in the turbulent stress tensor Tij we consider only the contribution of the 
turbulent kinetic and magnetic pressure as formulated in Eq.Q: 



(11) 



In a linear approximation (the weak mean-fleld case), Ki^2 = 1/3 (Kleeorin et al., 1996). In 
general, the turbulent stress tensor contains terms governing the differential rotation and 



the meridional circulation (Kitchatinov & Ruediger, 1995). We assume that these terms 



are stationary and do not contribute to the cross-helicity evolution. Therefore, we simplify 
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Eq.g as: 

d,{^) = -£:.(20 + W) + ^(B-V) fp^?+^ I -i^, (12) 

op \ Z/io / Tc 

Thus, the major sources of the turbulent cross-hehcity are due to the mean electromotive 
force, the mean vorticity, 2Q + W, and the gradients of the turbulent energy. For the 
simplest representation of the mean electromotive force, S ^ aB — tjtV x B + o 



L 

where the first term represents a-effect (here, a is pseudo-scalar) and tjt is the turbulent 
diffusion coefficient, we can write 



a,(u.b) = -(B.V)^u= + -| (13) 



- «(B-(20 + W))+/io77T(20 + W)-J-— , 

Tc 

The first two terms of this equation show that the mean cross-helicity is generated in the 
presence of the large-scale magnetic field that permeates the stratified turbulent medium. 
Another important contribution represented by the third term is due to the large-scale 
electric current and the mean vorticity. A similar equation was used by [Yoshizawa et al. 



(2000) for discussing a mean-field dynamo scenario based on the cross-helicity effects. 

A complete expression for the mean electromotive force was used by [Kuzanyan et al. 



(2007) who investigated the time-evolution of the cross-helicity produced by the mean-field 
dynamo in the solar convection zone. Their results showed that the cross-helicity sign 
alternates during the solar cycle with the amplitude of about ±2G kmS~^. The generated 
patterns of the cross-helicity depend on the structure of the mean electromotive force which 
is employed in the dynamo model. Their model, however does not include the stratification 



effects associated with the second term in Eq.(12). 



To study the distributions of the cross-helicity on time intervals much shorter than the 
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period of the solar cycle, we can neglect the time derivative in Eq.(13) and find: 



u-b ^^{B-V) l^pu^ + 1^ j - an (B . (2n + W)) + fioVrTc (2^ + W) • J (14) 

To estimate the magnitude of the cross-helicity in the near surface layers we use additional 

assumptions. First, we assume that the energy of the fiuctuating magnetic fields can be 

expressed through the kinetic energy of convective motions, = spu^^ and 6: = 1 is 

2/xo 

the energy equipartition condition. The convection in the surface layers is short lived in 
comparison to the period of the global rotation. For these condition the a-effect can be 
estimated as following a ^ tjtTc^ (6: + 1) V log ( pu'^ ) (see. 



Pipin 



2008). Then, the second 



term in Eq.(14) is much smaller than the first term because of <C 1 in the subsurface 



layer. Therefore we will neglect this contribution in the further study. For the further 
comparison we introduce the density stratification scale parameter, G = Vlogp, and 



rewrite Eq.(14) as follows 



u-b ^ r^T (^ + 1) { (B . G) + (B . V) log (u^) | + poVm (2^ + W) • J 



(15) 



If we neglect the effect of fiuctuating magnetic field, we find that the contribution of density 



stratification agrees with (Riidiger et al., 2010), but the contribution of the turbulent 



intensity stratification is larger by a factor of 2, in our case. The difference can be explained 



by the approximations made in calculation of Eq.(15). In particular, the contributions 



of (VU)^.^.^ (VB)^..^ are not included in Eq.(15). The main purpose of this approach 



is to demonstrate the basic physical contributions to the mean cross-helicity in the solar 



conditions. We have to point out that in the derivation of Eq.(15) we did not make 



assumptions that the turbulence inhomogeneity scale is much larger than the typical scale 



of turbulent fiows, £. This assumption is used in previous studies (Riidiger et al., 20101), and 
strictly speaking their results can be applied only to the case of a weakly stratified medium, 
1. 
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3.2. Calculation of u • b in the Fourier space 

Next, we estimate the contributions of terms in Eq.([2]) using the r-approximation 



made in the Fourier space (see, e.g., Rogachevskii & Kleeorin, 2000; Radler et al., 2003 and 



Brandenburg & Subramanian, 2005). BasicaUy this approach has the same shortcomings as 



Eq.(8). However, it aUows us to study the structure of the cross-hehcity in further detail. 



The derivations are explained in Appendix A. In addition, we calculate correlation UrK^ 
which can be estimated from the line-of-sight observations of velocity and magnetic field in 
a central part of the solar disk (cf., [Zhao et al. , 2011). 

Here we present the result of these calculations in the form that can be compared with 



Kleeorin et al. 


(2003 


) and 


Riidiger et al. 


(2010) 



u b 



Vt 



(2 + 3e) (B ■ G) 



(16) 



+ ^"^"^ ^! ^ (B • V) log (uA + rjTT, {e + 1) fXo • J) 



+ 



2 

VtTcHo 



(3 + s) J ■ W + ^ (23s + 5) (VU)^^^^.^(VB)^^_^.^, 



10 



where £ is a ratio of the kinetic and magnetic energies of the turbulent pulsations, 
rjT = u^Tc/S is turbulent diffusion is coefficient, and Tc is a typical convection turnover 
time. The structure of Eq.(16) corresponds to Eq.(|2]) though contrary to Eq.([l5|) the 
contributions of density stratification (G) and turbulent intensity stratification (Vlog ^u^ 
are decoupled. For the other coefficients we find = /Uo^^^ {3 + s), K4 = fXoi]TTc {e + 1) and 



AC5 = 2^ {23e + 5). Also, we see that the contribution of the stratification effects in Eq.(15) 



agrees but have slightly different coefficients i^i^sl ys (1 + e), and that the contributions 
due to the electric currents are in agreement. Note, that here similarly to Riidiger et al.| 



(2010) in derivation of Eq.(16) we used the scale separation assumption, £G <^ 1 (weakly 
stratified medium). 
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For the correlation of the radial components we obtain: 



JW 



6 



(17) 
(18) 



(5 - 2e) VrBr 



(Urbr) 



5 



4:er]TTc 



JW 



+ 



10 
70 



163£ + 41)V{i!7,}V{i%, 



(19) 



where, e^^^ is a unit vector in radial direction. Both Eq.(16) and Eq.(17) generalize the 



previous results of Kleeorin et al. (2003) and Riidiger et al. (2010) by including the effects 



of the large-scale electric current and velocity shear. In our models we find that the toroidal 
component of the large-scale magnetic field is much stronger than the poloidal component. 



We also discard the effect of meridional circulation in Eq.(19). Therefore we define the 
product of the magnetic and velocity strain by formula: 



where Vr^e are the covariant derivative components, B^j^ is the toroidal magnetic field and 
U = r ^iiiO {Vt (r, 6) — f^o) is the large-scale shear fiow due to the differential rotation. 



4. The Urhr patterns by dynamo models 



In the next section, we illustrate contributions of each term to the correlation of 



the radial components, Urhr^ given by Eq.(17) for two types of solar dynamo models: the 



mean-field dynamo distributed in the convection zone and the benchmark fiux-transport 



dynamo. The distributed dynamo model is described in detail by Pipin & Seehafer (2009) 



and Seehafer & Pipin ( 2009| ), hereafter, PS09 and SP09, respectively. The benchmark 
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model was constructed by Jouve et al. (2008). We consider the large-scale magnetic field 



produced by the models as an input parameter for estimation of Urbr in Eq.(17). The 
dynamo equations are given in Appendix B. 

The distributed dynamo model includes the meridional circulation with the geometry 



of fiow which is similar to that used by Bonanno et al. (2002). The maximum velocity of 



the meridional circulation is fixed to 10 ms ^. The turbulent generation effects include the 



a— effect (Parker, 1955) and the H x J effect (Radler, 1969). The internal parameters of 



the solar convection zone are given by Stix (2002). The integration domain is above the 



tachocline, from O.Tli?© to O.972i?0 in radius, and from the pole to the pole in latitude. We 



use the differential rotation profile given by Antia et al. (1998). The internal parameters 



of the model, including the turbulence parameters and parameters of stratification of solar 
convection zone, distribution of the angular velocity and the geometry of the meridional 
circulation are shown in Figure [!} 



For comparison we examine the UrK patterns of the benchmark fiux-transport dynamo 
model. This type of dynamo models is one of those widely discussed in the literature 



(e.g., Bonanno et al., 2002; Dikpati et al., 2004). We keep the same formulation of the 



benchmark dynamo model as described by Jouve et al. (2008), although we use the same 



type of rotation law as in the previous subsection. The integration domain in this model is 



from O.65i?0 to O.97i?0. The calculation of the UrW is carried out from O.71i?0 to O.97i?0. 
The magnetic field in this benchmark model is measured in the units of Bq. Therefore the 
calculated cross-helicity is scaled with the magnetic field strength Bq. The other parameters 
of the model are: the reference turbulent diffusivity, 5 • lO^^cm^s"^ ((7^ = 0.005), the 
magnetic Reynolds number, Rm = 700, the coefficient of the poloidal magnetic field source 



term, Cs = 30 (in notations of Jouve et al., 2008). 



We separate contributions to UrK into three parts: (urW) , {^rW)^j and (urbr) 
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and associate the first term of Eq.(17) as a contribution of the stratification effects, while 
the second and third terms represent contributions of the large-scale electric current. In 
the calculations we assume the equipartition between the kinetic and magnetic energies 
of fiuctuations, £ = 1. Modeling the cross-helicity we found that if we take the surface 



values for the model quantities in Eq.(17) then we obtain that the stratification effects are 
dominant. This is true for the both types of dynamo models. However, if we integrate 
the cross-helicity from the surface down to deeper layers, say down to O.9i?0, then the 
contributions of the large-scale current becomes comparable to the stratification effects 
and even greater. The dynamo butterfiy diagrams and the simulated patterns of the 
cross-helicity variations are shown in Figures [2] and |3j In the simulations, the antisymmetric 
modes of the toroidal magnetic field are dominant. Therefore we show the butterfiy diagram 
only for the northern hemisphere. The butterfiy diagram of the toroidal magnetic field 
is shown for the bottom of the convection zone (from where the magnetic field of active 
regions presumably erupts), while the radial magnetic field patterns are shown for the top 
of the integration domain. In the first model the toroidal magnetic field strength varies with 
the amplitude of 2.5kG at the bottom of the convection zone. The radial magnetic field has 
the maximum of the magnitude variations at the poles, where it varies with an amplitude of 
40G at the surface. The cross-helicity varies with amplitude ±10 G • km • s~^. In the second 
model, the magnetic field strength is measured in the units of 5o, which is assumed to be 
the maximum of the typical strength of the fiux tubes accumulated just beneath the solar 



convection zone (Dikpati & Charbonneau, 1999). The cross-helicity in the fiux-transport 
model is scaled with Bq as well. In this model, Bq should be about lO^G in order to obtain 
the cross-helicity values similar to the previous model. The main reason for such high Bq 
value is, of course, the low magnetic diffusivity used in the model. As shown by [DikpatT 



et al. (2004), the fiux-transport models can be tuned for high magnetic diffusivity in the 



near surface layers. In both cases, one can see that the part of the cross-helicity, which 
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results from the density and turbulence stratification effects ( panel (b) of the Figures 2 and 
3), follows the evolution of the radial magnetic fields while the large-scale current produces 
a more complicated pattern of the cross-helicity of comparable magnitude, panels (c) and 
(d). Both models give somewhat similar predictions for the cross-helicity patterns. The 
most obvious difference is that the second model predicts a very strong amplitude of the 
cross-helicity near the poles. This effect results from the strong polar magnetic field which 
varies in the range of ±0.3Bq while the toroidal magnetic field in this model varies in the 
range of ±3Bo. In Figure 3, in order to resolve the fine structure of the sunspot formation 
zone we restrict variations of B^. within the limits of ±0.01Bq. Comparing with observations, 
the polar magnetic field seems to be too strong in the benchmark fiux-transport model. 



Perhaps this model needs further tuning in a way similar to given by Dikpati et al. (2004). 



Looking at the contribution of the cross-helicity produced by the large-scale current we 
find that in the distributed dynamo model the surface cross-helicity pattern correlates with 
the variations of J^, and in the fiux-transport model both components of the current are 
important. The large-scale distribution of electric currents can be deduced from synoptic 



magnetogram (Pevtsov & Latushko, 2000; Wang & Zhang, 2009). For the current theory it 



would be interesting to compare the spatial-temporal patterns of the cross-helicity and the 
large-scale current. These patterns will evolve in phase if the contribution of the electric 



current in Eq(16) dominates 



If we consider the observational problem of finding the UrK^ we have to choose suitable 



spatial and temporal scales for averaging (in the spirit of Zhang et al., 2010), in order to 



distinguish between the processes of the "local" (noted in Mason et al., 2006; Matthaeus 



et al., 2008) and "global" alignments of velocity and magnetic fields. We can take as a 



working hypothesis that for larger scales most contribution to the observed cross-helicity 
comes from deeper layers of the solar convective zone. This assumption is consistent with 
our approximations but requires further investigations. In particular, our study neglects 
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turbulent fluxes of the cross-helicity. This assumes that the spatial scale for averaging 
exceeds the typical scale of the turbulent velocity variations. 



5. Discussion 

In this paper using the mean-fleld magnetohydrodynamics framework we derive an 
equation for the cross-helicity evolution in the presence of the large-scale magnetic flelds 



and flows. This equation (see Eq.(12)) can be used to construct the mean-fleld dynamo 



model based on the cross-helicity effects (see, Yoshizawa et al., 2000; Sur & Brandenburg 



2009). It is also can be used to estimate the cross-helicity distribution in subsurface layer 
of the Sun. The structure of the cross-helicity source terms in the evolution equation was 
studied on the base of the MHD equations governing turbulent fluctuations of velocity 
and magnetic flelds in stratifled media in the presence of the large-scale magnetic flelds 
and shear flow. We calculated the contributions of various effects to the cross-helicity 
pseudo-scalar. For comparison with solar observations we calculated the correlation of 
the radial components of the turbulent velocity and magnetic fleld and estimated these 
contributions for two different types of the solar dynamo. 

Our results show that the turbulent cross-helicity reflects internal properties of the 
dynamo processes and allow us to discriminate among different solar dynamo models. For 
analysis of the observational data involving the cross-helicity it is important to consider 
averaging over small and large scales. If cross-helicity data are averaged over a small portion 
of the solar surface and short time periods, the results may contain mainly the effect of 
the density stratiflcation and the vertical component of the large-scale magnetic fleld as 



shown by Riidiger et al. (2010). However, this type of averaging may meet some practical 
problems because of the difflculty to distinguish between the processes of the global and 
local alignment, of the turbulent velocity and magnetic flelds. For instance a similar kind 
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of alignment for velocity and vorticity is known as kinetic helicity (see Pelz et al., 1985). 



However, for different scales of convection motions near the surface level one may observe 



different signs of the kinetic helicity at different depths (Zhao, 2004). 



If the cross-helicity data are averaged over large spatial and temporal scales, then other 
effects related to the solar dynamo may become important (see Eq.(p!2|)). In particular, 
the cross-helicity patterns depend on the distribution of large-scale electric currents and 
velocity shear. In this case, the analysis becomes more complicated, and the large-scale 



transport of cross-helicity represented by T (see Eq.(lO)) has to be taken into account. 



Furthermore, the cross-helicity equations (16) and (|17[), may become invalid in the deeper 



layers of the Sun due to the strong Coriolis force effects [Vt^ ^ 1). 

Our study is done in the anelastic approximation and does not take into account the 



density fluctuations effects. However, it is known (e.g.,|Kichatinov fc Pipin, 1993) that in the 



turbulent coefficients the density fluctuations effects comes with factor gv? /C^ ^ v? /\Hp\^ 
and may become important. Therefore, we anticipate further development of the present 
theory regarding the turbulent transport and compressibility effects. In this paper, we 
have shown that the spatially-temporal patterns of cross-helicity depend on details of the 
dynamo mechanism. This suggests to use the cross-helicity tracer as a diagnostic tool for 
solar dynamo models. The forthcoming analysis of space and ground based observations 
towards direct determination of the cross-helicity challenges the future theoretical studies. 
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Fig. 1. — Internal parameters of the solar convection zone: a), the density scale Hp = 
(solid line), and the turbulent velocity scale Hu = (Vlogt^^)""^ (dashed line), the distance 
is measured in units of the solar radius; b) the typical turnover convective time, Tc (solid 
line), the coefficient of the turbulent diffusivity, rjT^ (dashed line), the RMS of convective 
velocity (dash-dot line); c) the contours of the constant angular velocity are plotted for 
the levels (0.75 — 1.05)1^0; d) the geometry of meridional circulation, the length of arrows 
are proportional to the circulation speed. 
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Fig. 2. — The time-latitude ("butterfly") diagrams for the distributed dynamo model. Con- 
tours in each plot the toroidal magnetic fleld at the bottom of the convection zone. The color 
background shows: a) the large-scale radial magnetic fleld at the surface; b) the cross-helicity 



component (urK) ; c) {urhr)^j ; d) {urhr) deflned in Eq. (17 18). 
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Fig. 3. — The same as in Figj2]for the flux-transport dynamo model. 
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Appendix A. Calculation of the cross-helicity in r-approximation in Fourier 

space 

Our approach for calculation of the cross-helicity is presented in detail by |Pipin| 



(2008), hereafter, P08. The important steps in procedure were described earlier (see, e.g., 
Rogachevskii & Kleeorin, 2000; Radler et al., 2003 and Brandenburg fc Subramanianl , 2005). 
The calculation is based on the equations governing the evolution of the fluctuating velocity 



and magnetic flelds in a rotating coordinate system (see, also, Rogachevskii & Kleeorin 



2000; 


Radler et al. , 


2003; 


Pipin, 


2008) 



dh 

dnii 



:Vx(u X B+U X h)+riV^h+&, 

2 



dt 



-2 (O X m) . -Vi\p' 



G • m) z/ + 



(b-B) 

2fi 



+ —Vj (Bjbi + Bibj) - Vj {Ujrui + Uirrij) 



(20) 



(21) 



They can be obtained from Eqs( |3[4[ ) by splitting velocity and magnetic flelds into the 
mean and fluctuating part (see Eq.Q). Here, 0, and ^ stand for unspecifled nonlinear 
contributions of the fluctuating flelds, m = pu, G = Vlogp is the density stratiflcation 
scale of the media, the fluctuating pressure component, ft is the angular velocity 
responsible for the Coriolis force, U is the mean flow velocity, weakly varying in space, and 
f is a random force driving the turbulence, r] and u are the microscopic magnetic diffusivity 
and viscosity, respectively. The further calculation is done as follows. First, it is convenient 



to write equations (|20p and (|21|) in the Fourier space: 

= u,/ 



%(z - q) ( (q) - mi{z - q) (^ ) (q) 



d<22) 



/ 



6,(z-q)V,(q)-6,(z-q)y,(q) 



dq + (J5j. 



dt 



+ uz^ + iu {Gz) ) mi{z) = fi + di - 2 {nz){zxni). 



(23) 
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ITT, 



mi{z - q)V/(q) + m/(z - q)Vz(q) 



+ -nif{z)zi I bi{z-q)Bf{(i) + bf{z-q)Bi{q) 
A* 



dq 

dq, 



where the turbulent pressure was excluded from (21) by convolution with tensor 



7rij(z) = 5ij — ZiZj, Sij is the Kronecker symbol and z is a unit wave vector. The equations 
for the second-order moments that make contributions to the cross-helicity tensor can 
be found directly from ( 22|23 ). As the preliminary step we write the equations for the 
second-order products of the fluctuating fields, and make the ensemble averaging of them, 

d 



dt 



rhi (z) bj (z') 



r/i^.(z, z') - {r]z''' + 1/^2 + (Gz)) (mi (z) bj {z'] 
rhi (z) rhjiz' - q)) (q)- 



(24) 



{rhi (z) "i/(z' - q)) 



Bi 



(q) 



dq - 2 (Oz) ennZi (rn„(z)6j(z') 



+ 



mi {z)bi{^' - q) ) Vj{q) - (mi (z) 6,-(z' - q) ) Vi{q) dq 



i7r,/(z)^; J [(^fhiiz - q)6, (z')) F/(q) + (m/(z - q)5, (z')) F;(q) 
k{z - q)5, (zO) Bf (q) + (bf{z - q)6, (z')\ Bi (q) 



d_ 
dt 



{rhi (z) nij (z')) = -2 (Oz) Sii^f; {rhn{z)rhj{z')) - 2 (Oz') £j7„2;; {rhi{z)rhn{z')) 



dq 

dq, 

(25) 



mif{z)zi {mi{z - q)mj {z'))Vf{q) + {fhfiz - qjrhj {z'))Vi{q) dq 



+ -'Kif{z)Zi 



at 



>,(z)5,(z') 



rhi (z) m;(z - q)) V f{(\) + {rhi (z) m/(z - q)) ^((q) 

bi{z - q)m, (z')) (q) + (6;(z - q)m, (z')) (g) 

rhi {■z)bi{z - q)^ % (q) + (rhi {z)bf(z - q)^ Bi (q) 
Th1j{z, z') - V {z'^ + z^ + i (Gz) + i (Gz')) {nii (z) rhj (z')) , 
r/^^/z,z')-(77^'^ + r7^2) (5.(z)6,(z') 



dq 
dq 
dq 
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+ 



k (z) m,{z' - q) ) ( y ) (q) - {k (z) m{z' " q)) (y ) (q) 



?fii(z - q)bj (z')) (q) - (mi{z - q)bj {z'] 

h{z)liz'-q))Vj{q)-(k{z) (z' - q) \ FKq 



P 

dq, 



(q) 



dq 
dq 



where, the terms T/i|^'^'^^ involve the third-order moments of fluctuating flelds and 
second-order moments of them with the forcing term. To proceed further we introduce the 
double Fourier transformation of an ensemble average of two fluctuating quantities, say / 
and taken at identical times and at the different positions x, x', is given by 



/(x)^(x') = 1 1 /(z)^(z')e^(^-''+^'-')d3zdV. 



(27) 



In the spirit of the general formalism of the two-scale approximation (Roberts & Soward 



1975) we introduce "fast" and "slow" variables. They are deflned by the relative r = x — x' 



and mean R = (x + x') /2, coordinates respectively. The fast and slow variables in the 
Fourier space correspond to two wave vectors: k = (z — z') /2 and K = z + z^ Then, we 



deflne a correlation function of f and g obtained from Eq.(27) by integration with respect 
to K, 



$ (/, ^, k, r) = I / (k + I) ^ (I - k) e^(^-^)d^ 



(28) 



For further convenience we define the foUowing second-order correlations of the momentum 
density, magnetic fluctuations, and cross-correlations of the momentum and magnetic 
fluctuations: 



Vij (k, R) 
p2u2 (R) 
hij (k, R) 
b2(R) 

Xij (k, R) 



$(mi,mj,k, R), 
Jvii (k,R) d^k, 
$(6i,6j,k,R), 

hii (k, R) d^k, 
$(mi,6j,k,R), 



(29) 
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We now return to equations (24), (25) and (26). As the first step, we approximate the 



Th\j'^''^^ -terms by the corresponding r relaxation terms of the second-order contributions. 



Th. 



ij 



Th> 



(v) 



Th 



(h) 

ij 



rhi (z) hj (z') 



{rhi (z) rhj (z')) - (m, (z) rhj (z')) 



(0) 



6,(z)6,(z'))-(6. (z)6,(z') 



(0) 



(30) 
(31) 

(32) 



where the superscript denotes the moments of the background turbulence. Here, 
Tcis independent on k and it is independent on the mean fields as well. As the next 
step we make the Taylor expansion with respect to the "slow" variables and take the 



Fourier transformation, (28), about them. The details of this procedure can be found in 



(Brandenburg & Subramanian, 2005). Furthermore, we consider the high Reynolds number 



limit (relevant for the solar turbulence) and discard the microscopic diffusion terms. Then 
we arrive to the evolution equation for the cross-helicity tensor (see, P08): 



dt 



-i(Bk)|^-M+l^ 



_|_ 'JUL 

p 



hi 
Po 



(Bk) dvii (GB) 1 . . _ . 

+ ^ ^-GsT^-^ ^-Vii kikfBfihu 

2p ' dks 2p ^0 



+ -GiViiBj + — '- 
P Po 



ViiBj^i k,j 



P 



kiBij d 
2 dkf 



P 



hj£ 
Po 



2k 



2kikf ( k • V 



kiVf - kfVi 



ki k i 



- 2 {Qkj kpEipikij - 2- \ Q.kj kpSipi ykV j kij 

i Q 
+ ^^ipi ( (^k^ VpXij + kp (1^ V) ki^ +kiU f^i 



(33) 
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where k is a unit wave vector, U f^i = V/[// and Bf^i = ViBf. The expressions for the 



Fourier image of cross-hehcity pseudo-scalar can be obtained from Eq.(33) if we assume 
quasi-stationarity (discarding the time-derivative terms). Then we solve the obtained 
equation in the linear approximation with respect to the large-scale fields and the Coriolis 
force and obtain the cross-helicity tensor kij. After this we calculate the cross-correlations: 



ub 



p J 



(34) 
(35) 



where e[^^ is the unit vector in radial direction. For the final results we have to define the 
spectra of the background turbulence. We consider the case of stationary and quasi-isotropic 



background turbulence(Kichatinov, 1987), 



p^E {k, R) 
SttP ■ 



(36) 



where 7r^j(k) = 6ij — kikj^ 6ij is the Kronecker symbol, the spectral function E{k^ R) defines 
the intensity of the velocity fiuctuations: 



^ / 



47rA;2 



(37) 



Similarly for the magnetic fluctuations of the turbulence, by the spectral function B{k, R) : 



b2 



(38) 



The anisotropy in Eq.(36) is induced by inhomogeneity of the turbulence. However, the 
scale of inhomogeneity is assumed to be small compared to the typical scale of turbulent 
fiows. Therefore, the turbulent anisotropy is small and the given model of the background 



turbulence is called as quasi-isotropic (Kichatinov, 1987). We omit the further cumbersome 



derivations, which are explained in detail in P08. 
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Appendix B. Dynamo model equations 

The evolution of the axisymmetric magnetic field {B being the azimuthal component 

of the magnetic field, A is proportional to the azimuthal component of the vector potential) 

is governed by the following equations: 

dA . , Uesine dA OA 

— = rsmOt^^ t/^— (39) 

ot r ofi or 

rsineS^ = (jl^^ Co^r^TGrnsmOB f['} (40) 

d^A 



+ 
+ 



(1 - ij?) d'^A 



dB _ rdndA _dndA\ _d{rUrB) 

dt \ dr dfj, 5/i dr ) dr 

sin 6 dUaB 1 drSa sin 9 dSr 



sin^5£'r sin^ 5 / , / ^(rf) 



r d/j, r dr r dji 
+ 2/<«.(l-,^)))^ (42) 

fj''^ sin 6^ ^ 9 , , . ^ 
r oil 

I drSe ^ 1^/.- ./ ^^(^) , c.Aa) 2\\9rB 
r dr r dr 



(*l*,(/r + 2/<°V0)^^+ (43) 



+ 



2/i sin ^ a ^ ^(a) , d 



1 a 

-r 



r 9r 

where /i = cos 6^, 6^ is co-latitude. We use the same notations as in P08. Functions 
fi^2^'^ depend on the Coriolis number fi* = 2tc^o^ is a typical correlation time of the 
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turbulent convection; functions cpf^ and il^r^^s describe magnetic quenching and depend on 
/? = B j \fji^p^ ^ parameters Cc^, control the contributions of the a and Q x J effects. In 
notations of P08, ipr] — 4^3 + 4^2 — 20i and ips — ^t"^ + 02^^- We introduce parameter to 
control the turbulent diffusion coefficient, t/t = Following to the hnear analysis 

of PS09 and SP09 we choose Cs = 0.02, = Cs/2, C^^ = 0.1. The dynamo equations with 
these input parameters give a solar-type dynamo solution. 



